Quasiparticles in the Kondo lattice model at partial fillings of the conduction band 
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I. INTRODUCTION 

The Kondo lattice model (KLM) has been a matter of 
constant interest for more than the last three decades. In 
two and three dimensions it is one of the common models 
to describe heavy-fermioni physics and is also a possi- 
ble candidate for high-T c superconductivity^ Our moti- 
vation to study the one-dimensional KLM2 is threefold. 
Firstly it has been shown^ that the spin polaron, which 
is a quasiparticle of the KLM, plays an important role 
in noncquilibrium transport in a quantum wire coupled 
to a ferromagnetic spin chain; our method provides the 
possibility to investigate the quasiparticles of the model. 
The spin polaron might also play an important role in the 
electron spin decay process^ in quantum dots induced by 
the hyperfine interaction due to nuclear spins. Secondly 
it might be helpful to understand the one-dimensional 
model in greater detail to assist investigations in higher 
dimensions. And lastly the model has become interest- 
ing for the description of mesoscopic systems, like carbon 
nanotubcs filled with fullcrcncs or cndohedral fullercncs, 
so called peapods£ The aim of this work is to expand 
on the understanding of the spectral properties of the Id 
KLM. We show, by means of the density-matrix renor- 
malization group&i (DMRG), that persistent quasiparti- 
cle states exist, namely the spin polaron states, and ex- 
trapolate their lifetimesand their weights. Furthermore 
we calculate the quasiparticle dispersion relation. For the 
case of half-filling we show that our results qualitatively 
agree with the results of a strong coupling expansion in 
Ref. H. We compare dispersion relations and confirm 
the existence of a critical coupling constant at which the 
effective quasiparticle mass diverges for large momenta. 

The KLM (Fig. QJ consists of a lattice with one local- 
ized f-electron on every of the L lattice sites, which do not 



interact with each other and a band of itinerant conduc- 
tion electrons of finite filling n, coupled to the localized 
electrons by an on-site Heiscnbcrg exchange interaction. 
The Hamilton operator of the Id KLM is given by 

L L 

H=-tJ2Y, (cU+i.+h-c) + J^Si-Si, (1) 

i=l cr=T,i i=l 

where t is the hopping parameter, c^ CT generates an elec- 
tron at site i with spin a and Si (sj) are the spin operators 
of the localized (conduction) electrons at site i, respec- 
tively. J is the Kondo coupling constant; we will consider 
only J > here, i.e. the antiferromagnetic coupling case. 
With k we denote the quasimomentum in the following. 

In principle, the Id KLM supports three phases, de- 
pending on the filling n and on the coupling J: A ferro- 
magnetic, a paramagnetic and, at half-filling (n = 1) 
only, a spin liquid phase, see Fig. [5] At half-filling 
of the conduction electron band the model is best un- 
derstood and early works using large- TV-methods^ and 
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FIG. 1: (Color online) The Kondo lattice model. The con- 
duction electrons are depicted in the upper row (green) and 
the localized electrons are depicted as bolt arrows in the lower 
row (red). 



the Gutzwiller approximatio n 10 ! 11 revealed that the mag- 
netic properties are due to the competition of the RKKY 
interaction and the formation of Kondo singlets, where 
such a singlet is a conduction electron forming a spin 
singlet with a localized electron. Due to half-filling, 
the electrons induce an effective RKKY interaction be- 
tween the localized spins, which forces antifcrromagnctic 
power-law correlations in the ground state. The oc- 
curence of RKKY oscillations or 2fcp-oscillations could 
be confirmed in Ref. [13 using DMRG. By means of ex- 
act diagonalizationi^ and quantum Monte Carloi^ it was 
shown that the ground state is spin- and charge-gapped 
and that it is a singlet of total spin. Therefore the ground 
state can be associated for all J with the universality 
class of spin liquids. There has been a controversial dis- 
cussion about the size of the Fermi volume (which is a 
single line in one dimension), whether it is small and 
therefore the Fermi wave vector is fcp = or whether 
it is large and therefore kp = ^ (n + 1). While a small 
Fermi volume would correspond to only conduction elec- 
trons contributing to the Fermi volume, a large Fermi 
volume would mean that the localized electrons also con- 
tribute to the Fermi volume. The idea of a large Fermi 
volume is borrowed from the periodic Anderson models 
There the f-electrons can move back to the conduction 
band and therefore contribute to the Fermi volume. The 
KLM can be derived from the periodic Anderson models 
in the case of large Coulomb interaction, where only one 
localized electron per site is allowed and other occupa- 
tions are fully supressed. This gives rise to the ques- 
tion whether the Fermi volume is also large in the KLM. 
Lately the authors of Ref. d could argue within a strong 
coupling expansion and from the evaluation of the con- 
duction electron density that the Fermi volume in the 
case of half-filling is small. In the same work, Ref. the 
quasiparticle dynamics of the half filled KLM have been 
examined as well. It has been possible to calculate the 
quasiparticle dispersion relation to good accuracy, where 
the quasiparticle mass has been found to diverge around 
k w 7r for J > J c 0.50 ± 0.02. Therefore the quasipar- 
ticles behave like nearly localized f-clcctrons due to the 
strong correlation of the conduction and localized elec- 
trons. 

In the limiting case of vanishing conduction electron 
density it could be rigorously showni£ by both ap- 
plying the Perron-Frobenius theorem and later exact 
diagonalizationii that the KLM is ferromagnetic for all 
J. Importantly, Sigrisl^ could show, that the quasiparti- 
cle of the Kondo lattice model is the spin polaron, which 
corresponds to an excited state separated from a con- 
tinuum of scattering states. Representatives of the con- 
stituent elements of the spin polaron are shown in Fig. 
[5J In Fig. [Ua) the localized spin lattice is completely 
ferromagnetic and the electron spin is oriented in the op- 
posite direction. Due to the antiferromagnetic exchange 
interaction the electron energy for quasimomcntum k is 
reduced. Via spin-flip processes, this state is coupled to 
the states shown in Fig. [U^b), where the electron spin and 



2 

(a) 

mi if tifiinii 




it mint nun 

q 

FIG. 2: ( Color online) Representatives of the constituent ele- 
ments of the spin polaron: a) Electron spin down; b) Electron 
spin up and one of the lattice spins down. 

one of the localized spins arc flipped and the momentum q 
has been transferred to the spin lattice. These are states 
of higher energy since the antiferromagnetic interaction 
can only reduce the energy if the two flipped spins are 
at the same site. The coupling leads to a level repulsion 
between the states of Figj2][a) and Figl^b), the energet- 
ically lower one corresponding to the spin polaron state 
and the higher ones forming the scattering states band. 
As proposed in Ref. d, the spin polaron state is expected 
to have a very long life-time if its energy lies outside the 
band of scattering states, so that it is protected against 
magnon absorption and emission processes. 

At partial band fillings n, ferromagnctism also survives 
in the strong-coupling limit where the KLM can be 
mapped to an effective Hciscnbcrg model with a ferro- 
magnetic exchange coupling. In this limit the formation 
of Kondo singlets, which move through the lattice, is suf- 
ficient to explain the occurence of ferromagnetism, but 
this does not exclude RKKY interaction, which might 
still play an important role. From exact diagonalization 
studies^ it follows that the KLM is ferromagnetic from 
a small but finite J(n) up to J — > oo for all n except 
half-filling. This raises the question, which mechanism 
drives ferromagnetism at intermediate couplings J and 
one proposal^ is that double exchange might be the cru- 
cial mechanism, where one conduction electron is respon- 
sible for screening several localized electrons. Screening 
is energetically effective in the antiferromagnetic KLM 
as long as J(n) surmounts a critical value and forces the 
localized spins to align in the same direction. 

At a certain J(n) the KLM approaches a second or- 
der transition^ by lowering J to a paramagnetic phase, 
where the spin polaroid might play an important role. 
The transition line has been calculated using exact 
diagonalizationi^ and has been refined later by means 
of bosonization^ The destruction of the ferromagnetic 
phase is described by a quantum random transverse- 



3 



field Ising Hamiltoniani^ Approaching the transition line 
from high J it has been propose d 20 i 23 that the large 
ferromagnetic cluster splits up in several small clusters 
each corresponding to one spin polaron. Just below the 
transition line the small clusters' direction of magnetiza- 
tion is not the same anymore for all clusters and leads 
to zero net magnetization. By means of DMRG the 
spin structure factor of the localized electrons could be 
calculated^ and it has been found that the size of the 
Fermi volume is small for very low J and becomes large 
approaching the transition line from lower J. From this 
one can conclude that near the transition line, the lo- 
calized electrons are incorporated in the Fermi volume 
and therefore spin polarons are formed. Lowering J the 
spin polarons are destroyed. The paramagnetic KLM has 
also been argued^i to belong to the class of Tomonaga- 
Luttinger liquids, 25 This is motivated by the gapless spin 
and charge excitation^ which also makes the model dif- 
ficult to handle with numerical methods using finite sys- 
tem sizes in this regime. From an analysis of Friedel 
oscillations, which are 2kp or Akp oscillations, the Lut- 
tingcr parameters could be determined for J > 1.8 and 
the Fermi volume has been found to be large. For very 
low J, RKKY or 2kp oscillations dominate the correla- 
tion functions of the KLM. This could be attributed 2 ^ 
to the backscattering of the conduction electrons at the 
localized electrons. 

In this paper we consider the spectral properties of the 
Kondo lattice model at partial band fillings. We will cal- 
culate the dispersion relation in the ferromagnetic phase 
for different Kondo couplings J and various fillings n and 
show that the spin polaron state exists. We are able to 
confirm the results of Ref. H at half-filling. In a second 
step we will show from the width of the spectral densi- 
ties that the lifetime of the spin polaron is very long and 
therefore the quasiparticle is persistent. We also examine 
the spectral densities in the paramagnetic phase and find 
that a spin polaron state can be excited, too. Its lifetime 
is smaller by several orders of magnitude than in the fer- 
romagnetic phase but the ratio depends very sensitively 
on the values of J, n and and the quasimomentum k. An 
interesting effect is found that the lifetime is maximal in 
the ferromagnetic phase if the quasimomentum is close 
to the Fermi volume. 

The paper is outlined as follows: In section [II] we will 
discuss the method, particularly how we calculate spec- 
tral densities, how we extract the lifetimes and how we 
extrapolate them. In section Mil we will present our re- 
sults. We will end up in a brief summary in section ITVl 



II. METHODS 

In this section we describe the methods used in our 
calculations. First of all we briefly discuss our DMRG 
algorithm. Secondly we describe the correction vector 
method, which we use to calculate Green's functions. At 
last we show how to calculate the lifetime of quasipar- 



ticles using the spectral functions we obtained from the 
Green's functions. 



A. DMRG 

For the calculation of ground states, we use a DMRG 
algorithm with abelian and non-abelian symmetries, 
whose implementation is based on a matrix-product for- 
mulation. We use open boundary conditions for all cal- 
culations. 



B. Correction vector method 

Applying the correction vecto r 25 ^ 1 i 28 i 29 method we ob- 
tain the spectral functions A(lu), where u is the energy. 
To calculate A(oj) : we need the retarded Green's function 
Ga{w + if}) = G\{uo + irj) + G^(u} + irj), where 



G + A {uj + irj) = (0)4* - 



1 



; A\0) (2) 



D + Eq + irj — H 

G^w + ir,) = (0\A — 1 -tt A^\0) (3) 

uj — Eq + irj + H 

are the two branches of the retarded Green's function and 
A is an arbitrary operator, |0) is the ground state with 
energy Eq and 77 > is an artificial broadening factor, 
which is needed to lower the lifetime of the excitation to 
avoid boundary effects due to the finite system size. The 
basic rule is to choose r\ > -j-, where c is the velocity of 
the excitation, but the minimal 77 is strongly depending 
on the model. 

In principle one would need to compute both branches 
of the Green's function to obtain the complete spectral 
properties. For the determination of life times the spec- 
tral weight of the quasiparticle is nearly completely con- 
centrated in one of the branches. Therefore we can ne- 
glect the other branch in this case. From now on, we 
will base all our arguments concerning the Green's func- 
tion on the +- branch. Concerning the spectral density 
the calculations for the — branch can be done similarly 
except for a minus sign. 

The correction vector is defined as 



|c(w + irj)} 



1 



uj + Eq + irj — H 



A\Q) 



and hence 



(w + Eq + irj - H) \c{uj + irj)) = A\0), 



(4) 



(5) 



where the groundstate |0) ist obtained from the preceding 
DMRG calculation. This leads to a non-hermitian system 
of linear equations for real and imaginary parts, which 
can be solved using the GMRES method^. The out- 
come is \c{lo + ir/)}, which allows to calculate the Green's 
function as 



G A (u + irj) = (0\A\c{uj + ir-/)) . 



(6) 



4 



The spectral density can then be obtained by applying 
the standard formula 



A(lo + irf) = Im Ga(lo + iff), 

7T 



where lo is assumed to be real. 



(T) 



with eo(k) = —2tcoska, where a is the lattice spacing. 
The one electron Green's function is then defined as 



G ka (uo + irj) 



1 



\Cka 



to + Eq + irj — H 



4J0). (io) 



C. Quasiparticle lifetime 



The self-energy T, a (k 7 Lo) is implicitly defined for the 
interacting system H as 



For the calculation of quasiparticle lifetimes we will 
limit ourselves to electronic systems. It is useful to trans- 
form the Hamiltonian into the fourier space according to 



Hence we obtain 



fc <T=T,1. 



(8) 



ff = EE (e (k)clc ka ) +jJ2Sk-S-k (9) 



G ka (uo + irj) 



1 



lo + irj - (e (k) - fi + T, a (k,Lo + irf)) '' 



with fj, the chemical potential. Note that rj appears also 
in the self-energy. This is nescessary, because lim^^o 
will not be carried out in the numerical calculations. In 
general, the self-energy is a complex function Ti a (k, lo) = 
R a (k,to) + iI a (k,Lo). Separation of real and imaginary 
part leads to 



G ka (u + irj) = 



1 



lo - (eo(fe) - M + Ra(k, to + irf)) + i(r/ - I a (k,Lo + irf}) 
I 



(11) 



We now assume, that the self-energy is continuous and 
only weakly depending on lo in the vicinity of a reso- 
nance lo 1(t = e (k) - [i + R a (k,Lo)\ iAj=LJ ^ : where Lo la is 
one out of several resonances, which are well separated 
to provide the correct determination of the lifetime of 
the quasiparticles (see the end of this section for the 
explicit extrapolation scheme). In addition we assume 
\I a (k,Lo)\ -c |eo(fc) — /i + Ra(k, lo)\ near to the resonance 
we are interested in, i.e. we assume long lifetimes, be- 
cause we are interested in these. This leads to 



I a (k,LJ + i V )^I^(k), 



(12) 



in the vicinity of the ith resonance. For the real part 
of the self-energy we apply a Taylor expansion at the 
resonance We find 

lo - (eo(fc) - H + Ra(k,LO + irj)) sa 



{LO - LO ia ) 1 



dR a (k 1 to + irf) , 



dRaik, lo + in) 
irj- 



dto 



dto 



\uj-\-irj—uii 



and define 



Otirr = 11 — 



dR a (k, lo) 



du> 



(13) 



Substituting this to Eq. (fTTj) the Green's function in the 



vicinity of resonance LOi a is given by 

i 

G k(J (lo + irj) = a ia - 



ia + i (v + a ia 
and the spectral function obtains the form 
A k<J (uo + irj) = 

r«(fc) 



i ( J\k) 



n (LO-LO i<T ) 2 +(r, + a i(T \l^(k)\) 



(14) 



(15) 



which corresponds to a sum of Lorentz distributions at 
the resonances uii a with a broadening of 



B^( V )= V + a i(7 jW(fc) 



(16) 



Hence the broadening computed with the correction- 
vector method is essentially the sum of the natural broad- 



ening a la 



li l \k) 



and the artificially introduced broad- 
ening rj and therefore B(r/) linearly depends on r\. Note 
that from the Lehmann representation of the spectral 
density one can find that I^(k) < 0. 

The broadened spectral density A ka (lo + irj) is a con- 
volution of the non-broadened spectral density A ka (uj) 
with a Lorentzian of width ?/: 



1 



A k(7 (uo + irf) = - / dLo'A ka (Lo) ■ — 

TT J (LO — lo'Y + rf 



(17) 
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FIG. 3: (Color online) The linear fit of B(rf) vs r\ is shown 
for three different quasimomenta, for a system with 48 sites, 
n = 0.125 and Kondo coupling J = 1. The data has been 
offset by r\. 
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FIG. 4: (Color online) Dispersion relations of the half- filled 
KLM. The comparison to 11th order perturbation theory is 
taken from The line for tj J = 0.8 is meant only as a guide 
to the eye. 



We now assume, that the spectral density consists of a 
sum of Lorcntz distributions, which are separated by non- 
lorentzian regions. The outcome of the convolution of 
two Lorentzians again is a Lorentzian, where the broad- 
enings behave additively. As the broadening corresponds 
to an inverse lifetime, we can define the lifetime of the 
quasiparticlc corresponding to the ith resonance as 



r = lim 



1 



V + a i0 
1 



(18) 



The limitations of this method are obvious. First of all 
the excitation must cause a lorentzian shaped peak in 
the spectral density. To be able to extract the broaden- 
ing of such a peak all other peaks must be energetically 
separated from this one. Thus we have to check, whether 
the conditions of our theory are fulfilled or not. We can 
check whether the spectral density has a lorentzian shape 
(see Fig. |5J) and B(rj) has to depend linearly on rj (see 
Fig.©. 

We have now two possible estimates for the lifetime of 
a quasiparticlc: 

1. Lower estimate: Use the inverse broadening B(r])~ 



directly (without 77 - 
smaller than B(Q)^ 1 , 
lower estimate. 



0). With (0 B^)- 1 is 
this is therefore a reliable 



2. 



Extrapolation: Calculate the broadening for sev- 
eral different 77. From this one can extract the lin- 
ear dependence of B(rf) on r] and £>(0) -1 gives the 
extrapolated lifetime. See also Fig. [3J 

It turns out that due to long lifetimes only extrapolated 
lifetimes are meaningful. Therefore in the next section 



we will discuss the results obtained by the second method 
only. 

From the Lorentzian fit of a single resonance peak of 
the spectral function in Eq. (|15[) one can also estimate 
the spectral weight of the corresponding excitation. 



III. RESULTS 




FIG. 5: (Color online) Dispersion of a KLM with 48 Sites, 
n = 0.125 and J = 1. The broadening is 77 = 0.2. The lower 
band is the bound quasiparticles band, while the upper one 
is the scattering states band. 



In this section we will present the results obtained us- 
ing the methods we discussed in the last section. First we 
will show the calculated dispersion relations considering 
a Kondo lattice model at half-filling and at partial fill- 
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ing. Our half-filling results show a qualitative agreement 
with the results in Ref. [||. In the second part we show 
several spectral functions and the calculated lifetimes, 
which leads to the conclusion that we find a bound po- 
laron state. The hopping parameter t is set to t = 1 in 
all calculations. 



A. Dispersion 

The half filled KLM serves as the touchstone of our 
method, where we can compare our results to those of 
Trebst et ali£, who did a strong coupling expansion up to 
11th order in t/J. We calculated the dispersion relation 
for different values of t/J, see Fig. 2]and used lattice sizes 
of 32 and 48 sites. The calculations have converged in 
the sense that we could not find any deviations between 
calculations of different system sizes. Our results show 
very well agreement to the results in Ref. [H for tj J = 0.4 
with small deviations for small k. By strong coupling 
expansion it is found that the band flattens out for k — > 7r 
around t/J c f» 0.50 ± 0.02 and therefore the effective 
quasiparticle mass diverges. This is also found by DMRG 
for a higher value of t/J c ~ 0.576 ± 0.002. As one can 
see, in Fig. [4j the strong coupling expansion dispersion 
relation has a minimum at k < ir for t/J = 0.6, which 
is not found by DMRG. DMRG exhibits this minimum 
again only for higher values of t/J. In Fig. 2] it is for 
instance shown for t/J = 0.8. Summarizing, the DMRG 
results match very well to the strong coupling expansion 
for low t/J but the agreement becomes worse for t/J> 
0.6. DMRG is the more reliable method in that regime, 
because it is non-perturbative and the error can be easily 
controlled by very small DMRG truncation errors. In this 
case it is easy to keep the truncation error reasonably low. 
We can confirm the physical picture established by Trebst 
et al., namely that the quasiparticlcs gain an enormously 
high mass, which is due to a growing correlation between 
the conduction and the f-electrons. The quasiparticles 
with high momenta therefore behave like almost localized 
f-electrons. 

Now we consider the dispersion relation of the KLM 
for partial band filling, sec Fig. [5] For the ground 
state calculation of the KLM with 48 sites, a filling of 
n = 0.125 and J = 1 we used about 100 DMRG ansatz 
states. The calculation of the correction vectors needed 
800 DMRG states to reach good convergence. In Fig. [5] 
and all other figures of dispersion relations we neglected 
the chemical potential, which would shift the lower band 
edge to nearly zero. One can distinguish two different 
bands. The higher one behaves like 2tcosk + 2t and 
can therefore be attributed to free electrons, which do 
not form bound states with the localized spins. From 
now on, this band will be referred to as the scattering 
states band in the sense that these excited states rapidly 
decay. The lower one of the two bands represents the 
states of the system which are formed by the conduction 
electrons bound to the localized spins, from now on re- 
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FIG. 6: ( Color online) Dispersion relation of a KLM with n = 
0.125 and J = 3.5. The lower band is the bound quasiparticles 
band, while the upper one is the scattering states band. 



ferrcd to as quasiparticle or Spin-polaron band. Contrary 
to the scattering states band this band consists only of 
one state, which is separated from the continuum (for 
large N) of scattering states and has a Lorentzian shape 
from which the lifetime can be extracted, which is very 
long in most of the cases, see Sec. IIII B[ Keeping the 
same filling n = 0.125, but raising the Kondo coupling 
constant J to 3.5, sec Fig. [HI the quasiparticle band 
becomes more separated from the scattering states, be- 
cause the quasiparticle state is now energetically lowered. 
This can be understood by a simple physical picture. For 
that we rewrite the exchange coupling of the Hamilto- 
nian as J2i [J Z S? sf + 4^ (S^s~ + S^s^ - )] and we now 
set J± = 0. The quasiparticle state of the KLM almost 
only consist of an electron with spin antiparallel to the 
localized spins. With respect to our modified exchange 
coupling this results in a lowered energy of J z /4 per 
electron. The scattering states also contain important 
contributions with an electron spin oriented antiparallel 
to the localized spins. This leads to a higher energy of 
J z /4. Therefore the energy difference between quasipar- 
ticle and scattering states is J z /2 and scales with J z . 
Taking also a finite Jj_ into account the quasiparticle en- 
ergy is even lowered more due to spin flip processes. The 
scattering states band has a similar shape as before, as 
expected. The weight of the quasiparticle states band is 
also increasing with J. This is also expected, because the 
state becomes energetically more favorable with increas- 
ing J. 

In Fig. [7] we show the dispersion relation of the quasi- 
particle of a system with J = 3.5 for three different fill- 
ings, n = 0.125, n = 0.25 and n = 0.75. The ground state 
is ferromagnetic in all cases. We conclude that even in 
the presence of many electrons the spinpolaron state can 
be clearly identified. 
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FIG. 7: (Color online) Dispersion relations of a KLM with 
J = 3.5 and three different fillings, n = 0.f25, n — 0.25 
and n — 0.75. The inset shows the second derivative of the 
dispersion relation, which is oc -nr, where m* is the effective 
mass. 
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FIG. 9: (Color online) Simplified phase diagram of the ID 
Kondo lattice model taken from [23|. The points mark the 
parameters at which extrapolated lifetimes have been calcu- 
lated. The lifetimes for points 1,2 and 3 and quasimomentum 
k = j§ are given directly in the picture by Ti, T2 and T3, the 
lifetimes for the other points are listed in tab. |T] 



B. Life time estimations from spectral functions 



In a further step we take a look at single spectral densi- 
ties for fixed quasimomentum fc, which provides the pos- 
sibility to calculate quasiparticle lifetimes of the bound 
quasiparticles and prove the existence of bound polaron 
states. We consider only the calculation of the extra- 
polated lifetimes, as described in section Hi CI and whose 
extrapolation scheme is shown in Fig. [3] Calculating 
extrapolated lifetimes this way we have to be very careful 
due to the assumption that the spectral density complies 
a lorentzian distribution. For the spectral density being 
a lorentzian the imaginary part of the self-energy has to 
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FIG. 8: (Color online) Spectral function of a KLM with 48 
Sites, n — 0.125, J = 1, k = ^ and r\ = 0.02. L, 1UJ0 (lu) is a 
Lorentzian fitted to the DMRG values. 



be very small compared to the energy of the resonance 
and it should not vary too much in the vicinity of the 
resonance. The expansion of the self energy leads then 
to a lorentzian function. Thus the spectral density is not 
lorentzian shaped in higher orders of the expansion and 
it has to be checked, see for example Fig. [8j whether it 
is good enough. Fig. [8J shows a spectral density for a 
KLM with 48 Sites, n = 0.125, J = 1 and quasimomen- 
tum k = jjj. The artificial broadening is set to 77 = 0.02. 
The number of data points obtained provides the possi- 
bility of a very precise fit of the lorcntz distribution. Fig. 
[5] sharply supports the assumption made in sec. Ill CI 
that the spectral density has a lorentzian shape, which is 
nescessary to calculate quasiparticle lifetimes. 

The extrapolated lifetimes are summarized in Fig. [5] 
and Tab. fl] There we can see, that the lifetime strongly 
depends on the parameters filling n and Kondo coupling 
constant J as well as on the quasimomentum k. The life- 
times in the ferromagnetic phase (this concerns the [n, J]- 
pairs {[0.125, 0.5], [0.125, 1], [0.25, 1], [0.29, 1], [0.75,3.5]} 
) decrease by approaching the paramagnetic phase by ei- 
ther lowering J or increasing n. For fixed and low quasi- 
momentum k it seems that the life-time decreases by in- 
creasing n (even if J is increased at the same time so 
that the distance to the paramagnetic phase is still large, 
compare e.g. the pairs [0.125, 1] and [0.75,3.5]). This in- 
dicates that the life-time is influenced by the presence of 
other quasiparticles, probably by an effective interaction 
between the quasiparticles mediated via the coupling to 
the localized spins. This is further substantiated by the 
dependence of the quasiparticle life-time on the quasimo- 
mentum k in the ferromagnetic phase. For k approaching 
the Fermi level, the lifetime increases, which is consistent 
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with the fact that the phase space for quasiparticle in- 
teraction becomes smaller close to the Fermi level. In 
contrast, clcctron-magnon interaction is expected to be 
more important for quasiparticles close to the Fermi-level 
because the energy of the spinpolaron comes closer to 
the scattering band. This effect can be seen in the para- 
magnetic phase for the pairs [0.375,0.5] and [0.75,2.5], 
where the lifetime decreases with increasing quasimo- 
mentum. Thus, in the paramagnetic phase, we conclude 
that electron-magnon interaction limits the life-time of 
the spinpolaron. Deep inside the paramagnetic phase 
at [n = 0.375, J = 0.5] the lifetime is short for all de- 
termined values of k. Therefore, as predicted earlier in 

Refs. mm there exist no persistent quasiparticles in 

this regime. 

We also extracted the spectral weight of the spin po- 
laron excitation from the Lorentzian fit and summarized 
them in Tab. [J in the second row of the respective k 
value. Considering the three numbered points of Fig. [S] 
we calculated the spectral weights 

1. 0.818004 ±0.000001 

2. 0.87119 ±0.00003 

3. 0.588 ±0.001, 

which do fulfill the expectation that the spectral weight 
should be significantly lower in the paramagnetic phase. 
The calculated weights are independent of 77 within the 
crrorbounds. They show a strong dependence on the 
quasimomentum (decreasing for growing k) in the fer- 
romagnetic as well as in the paramagnetic phase. This 
is expected because the spinpolaron states with higher 
value of k have higher energy and come closer to the scat- 
tering states. However, it is unexpected that the spectral 
weight is large for [0.375,0.5] and this still has to be ex- 
plained. 



IV. SUMMARY 

We have studied the one-dimensional Kondo lattice 
model at half-filling and at partial band fillings for vari- 
ous Kondo couplings J. At half-filling we could verify the 
results of Ref. [||. This includes the dispersion relation 
and the divergence in the effective quasiparticle mass. At 
partial band fillings we were able to show that in the case 
of ferromagnetism long living quasiparticle states exist, 
the spin polaron quasiparticle. The lifetime exceeds the 
lifetime of quasiparticle excitations deep inside the para- 
magnetic phase by several orders of magnitude. From the 
dependence on the quasimomentum we conclude that the 
dominant decay process is the spinpolaron-spinpolaron 
interaction in the ferromagnetic phase, and the interac- 
tion between electrons and spin waves in the paramag- 
netic phase. The weight of the spinpolaron state is very 



close to one even for special points in the paramagnetic 
phase. The results motivate the speculation that spin 
coherence can be significantly enhanced by coupling of 

TABLE I: Summarization of extrapolated lifetimes (first row 
of the respective k value) and spectral weight of the spin po- 
laron excitation (second row of the respective k value) for 
different fillings of the system, different coupling constants 
and different quasimomenta. The k values in paranthesis cor- 
respond to [0.75, 2.5]. 

^> 1 [0.25,1] [0.29,1] [0.375,0.5] [0.75,3.5] [0.75,2.5] 



(■|f) 831 ± 16 183 ± 3 67 ±9 
0.8211 ± 0.8292 ± 0.951 
0.0004 0.0008 0.006 



38.6 ± 0.8 25.7 ±1.6 
± 0.714 ± 0.674 ± 
0.008 0.005 



2540 ± 84 



14.7 ± 1.5 
0.82 ± 0.01 



2850 ± 
600 

0.6244 ± 0.6375 ± 
0.0001 0.0005 



36tt / 24tt 
40 V 33 I 



O(l0 3 ) 5.29±0.14 

0.5580 ± 0.3391 ± 
0.0008 0.0005 

O (10 3 ) 12.76 ± 
0.62 



0.3391 
0.0005 



± 0.190 
0.006 



electrons to magnons in the ferromagnetic phase of the 
localized spins. Even if the localized spins have a negligi- 
ble exchange interaction, the itinerant electrons can me- 
diate an effective exchange interaction which drives the 
local spins into a ferromagnetic phase at low enough tem- 
peratures (the critical temperature has been estimated 
recently and is within experimental reach. 31 i 32 i 33 ) As we 
have seen in this paper, the spin excitations in the fer- 
romagnetic phase can in turn form spinpolaron bound 
states with the itinerant electrons, increasing their life- 
time considerably. This effect persists in the presence of 
many electrons and becomes more efficient for quasimo- 
menta close to the Fermi level. It is an interesting ques- 
tion for future research to investigate the consequences 
for the spin relaxation and dephasing rates in the Kondo 
lattice model by studying the spin-spin correlation func- 
tions. 
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